Spin drag Hall effect in a rotating Bose mixture 
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We show that in a rotating two-component Bose mixture, the spin drag between the two different 
spin species shows a Hall effect. This spin drag Hall effect can be observed experimentally by 
studying the out-of-phase dipole mode of the mixture. We determine the damping of this mode due 
to spin drag as a function of temperature. We find that due to Bose stimulation there is a strong 
enhancement of the damping for temperatures close to the critical temperature for Bose-Einstein 
condensation. 
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Introduction — Electronic transport is one of the main 
topics of interest in condensed-matter physics, and an 
especially important phenomenon in electronic transport 
is the Hall effect. It was discovered already in the late 
nineteenth century by Hall. He observed that if a mag- 
netic field B is applied perpendicular to the current den- 
sity j through a conductor, the Lorentz force leads to a 
voltage drop in the direction perpendicular to both the 
current density and the magnetic field. This voltage is 
proportional to j x B, with a proportionality constant 
that depends only on the density of electrons and not on 
any other material parameters [1| . While the discovery of 
the Hall effect predates that of the electron, it is impor- 
tant for our purposes to note that electronic transport 
is in fact fermionic transport since it is mediated by the 
movement of electrons. By now many variations of the 
Hall effect have been found: the (integer and fractional) 
quantum Hall effects in which the Hall voltage is quan- 
tized @-Q , the spin Hall effect [3, @ , and the quantum 
spin Hall effect [7[. Recently, a spin Hall drag effect has 
also been proposed [8|. The latter is, as all spin Hall ef- 
fects are, due to spin-orbit interactions that play no role 
in atomic Bose gases if they are not externally introduced 
by applying an appropriate laser field configuration |9|. 
The proposal of Ref. [SJ is therefore physically very dif- 
ferent from what we discuss below. 

An important field of physics which connects few-body 
atomic physics with many-body and condensed-matter 
physics is that of cold atoms. The Bose-Einstein conden- 
sation of bosons at very low temperatures was already 
predicted by Einstein in 1924, but only observed directly 
in 1995 [10|. Over the past fifteen years, techniques have 
been getting steadily more refined, and it is now possible 
to make all sorts of degenerate atomic mixtures consist- 
ing of several spin states or of several different atomic 
species of either fermions or bosons. These mixtures are 
always trapped in optical and/or magnetic potentials, in 
which the atoms can be set into rotation by stirring with 
a so-called laser spoon |ll| . 

It is tempting to combine the above two fields to also 
get more insight into the physics of bosonic transport. 



At first sight this seems less than straightforward, since 
it is not possible to simply attach leads to a cloud of 
cold atoms to set up a steady-state transport current of 
atoms through the mixture. Moreover, in the cold-atom 
situation there are no obvious mechanisms that relax the 
particle current and give nonzero resistivities. However, 
we can use the phenomenon of spin drag as a bridge be- 
tween these two worlds 12Ml4l| . Spin drag was first pro- 



posed by D'Amico and Vignale by making an analogy 
with Coulomb drag between to electron layers 15|. It 
was later observed by Weber et al. [l6|. Whereas in the 
classic Coulomb drag experiment electrons are differenti- 
ated by the layer they occupy, in spin drag the spin of the 
electron is the relevant degree of freedom, i.e., electrons 
of one spin species drag along electrons of the other spin 
species. The resistivity created by this spin drag, which 
is a resistivity to spin but not to charge currents, typi- 
cally goes as pd oc T^ in electronic systems, with T the 
temperature. This is the trademark of a Fermi-liquid like 
behavior. 

In our earlier work |13| . we investigated the situa- 
tion in which the particles involved in the spin drag are 
bosons instead of fermions. We proposed an idealized 
set-up in which spin-1 bosons in the state {mp = +1) are 
accelerated along a torroidal trap by a time-dependent 
magnetic-field texture that creates a fictitious electric 
field. The atoms in state \mp ~ 0), which arc also 
present, do not feel this force but experience spin drag 
due to collisions with the other species. We found that 
due to the Bose enhancement of interatomic scattering, 
the drag resisitivity increases at lower temperatures. For 
the one-dimensional set-up considered, it in fact behaves 
as Pd oc T^^/^ for low temperatures, in strong contrast 
with the usual quadratic Fermi-liquid result. 

In this Letter, we discuss spin drag in a rotating Bose 
mixture. We consider the realistic situation of a three- 
dimensional Bose mixture with two spin components, 
present in equal numbers, just above the temperature for 
Bose-Einstein condensation. We first consider the homo- 
geneous case and look in linear response for steady-state 
solutions of the appropriate Boltzmann equation with a 



nonzero spin current. We find that the drag resistivity 
now becomes a 3 x 3 matrix with nonzero off-diagonal ele- 
ments that are proportional to the rotation speed, which 
represents a Hall effect. Indeed, these off-diagonal ele- 
ments are analogous to those found in the classical Hall 
effect and, in particular, do not depend on the specific 
collisional details of the mixture that determine the di- 
agonal resistivities, but only on the atomic density and 
external rotation frequency. 

As mentioned previously, such steady-state solutions 
no longer exist in the realistic situation that the atomic 
mixture is trapped in an external harmonic potential. 
The spin drag Hall effect can nevertheless be observed 
in that case by considering the collective modes. In par- 
ticular, we consider the dipole mode in which the two 
spin components oscillate out of phase with each other, 
because this mode obtains an orthogonal, i.e., a trans- 
verse component due to the spin drag Hall effect. More- 
over, the longitudional spin drag leads to damping of this 
mode, which makes it interesting to find out how the re- 
laxation rate of these modes depends on temperature. 
To obtain this, we again solve the Boltzmann equation 
for this specific case in linear response, and find that the 
relaxation rate shows a substantial increase as the tem- 
perature gets closer to the critical temperature. In three 
dimensions the relaxation rate does, however, remain fi- 
nite at the transition temperature. 

Spin drag Hall effect — To illustrate the spin drag Hall 
effect we consider first a homogeneous three-dimensional 
Bose mixture of two spin states, which we label |0) and 
|1), in the normal state. We assume that the bosons in 
spin state |1) couple to an external force F, and that the 
other spin state does not couple to this external force. 
(Generalizations to more than two spin species and dif- 
ferent forces are straightforward.) This would for exam- 
ple be the case if the external force is due to the Zeeman 
effect in a magnetic field, and if the two spin states cor- 
respond to the mp = and mp = 1 projections of an 
F = 1 hyperfine state. We assume that the system is 
rotating, which gives rise to a Coriolis force that is the 
equivalent of the the Lorentz force from the electronic 
Hall effect. 

The appropriate Boltzmann equation is 
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Here, /i (k, t) is the distribution function for the bosons in 
state |1). The Boltzmann equation for /o(k, t) is found 
by replacing /i <->■ /o and setting F to zero. Further- 
more, H = ilz is the rotation vector, and 17 x k gives the 
Coriolis force. The colhsion term, Tcoiiifo, fi], describes 
collisions of atoms with different spin. We will give its 
precise definition later on. There are of course also col- 
lisions between atoms with an identical spin but they do 
not play a role for the spin drag. 



We solve the Boltzmann equation by using the ansatz 
/i (k, t) = Nb (ek-mvi(t)/fi) ; with a similar expression for 
/o(k,t) in terms of Vo(i). Here, Nsit) = [e^^'-f") - l]-^ 
is the Bose-Einstein distribution function with /3 ~ 
{kBT)~^, as usual, the inverse thermal energy, fcs Boltz- 
mann's constant, and T the temperature. The single- 
particle dispersion is ek = h^k.^/2m with m the particle 
mass. The chemical potential /x is determined by the 
condition that the density of atoms is constant. The 
Boltzmann equation leads to the following equations of 
motion for the drift velocities Vq (t) and Vi (i) 



nm—— = Znmil x vq — 1 (vo — vi); 
dvi 



(2) 



nm = nF + 2nmri x vi + Ffvo — Vi). (3) 

dt 

Here, n is the particle density per spin. We assume this 
density to be equal for the two spin species. Again, note 
that generalizations of the above to spin and mass im- 
balanced systems are straightforward. 

In general the frictional spin drag is determined by 
the full nonlinear (vector- valued) function r(vo — vi) = 
/ dkhkYcowiNB (ek-mvo/n) > ^b (ek-mvi/fi)]/(27r)3. In 
the linear-response regime where the velocities are small, 
we make use of the fact that it can be approximated by 
r(vo — vi) ~ r'(0)(vo — vi). Note that here we make 
use of the isotropy of the collision term in the Boltzmann 
equation. We now introduce j = n(vi — Vq) the relative 
particle current, which up to dimensionful prefactors is 
equal to the spin current, and solve the above equations 
of motion for the steady state, i.e., dj/dt = 0. We then 
find that j = <T-F = p~^-F, which defines the con- 
ductivity and resistivity tensors a and p, respectively. 
Note that these are 3x3 matrices since the force and 
current are three-dimensional vectors. We find that the 
longitudinal resistivities pxx = Pyy = Pzz = 2r'(0)/n^, 
which are related to the spin drag relaxation time r via 
a Drude-like formula as pxx = m/riT. This relaxation 
time is the time scale on which the spin current decays 
due to collisions of atoms in different spin states. The 
Hall resistivities are given by pxy = —Pyx = 2mU,/n. 
All other components of the resistivity and conductivity 
tensors are zero. Like the Hall resistivity in electronic 
systems, the transverse components of the resistivity do 
not depend on the specifics of the processes that lead to 
a nonzero longitudinal resistivity, but only on the den- 
sity and strength of the Coriolis force. The longitudinal 
resistivity, however, that determines the dissipation of 
the relative momentum current via frictional spin drag, 
depends on the inter-spin-species collisions. 

Collective modes — To implement the spin drag Hall 
effect in a realistic cold- atom experiment, we have to take 
into account the effects of the trapping potential. Steady- 
state current are now no longer possible. In this situ- 
ation the collective-mode spectrum of the mixture pro- 
vides an experimental method to determine the spin drag 



resistivities. We consider a harmonic trapping potential 
Vtrap(x) = ma;o(x^ + y'^)/2 + mujlz^ /2 with radial trap- 
ping frequency wq and axial frequency w^. We now have 
for the Boltzmann equation 

where ^(x) = Vtrap(x) — mVi^{x'^ + j/^)/2 includes the 
centrifugal force. The equation for /g is again found by 
replacing fi ^ fa- 

We solve this inhoniogcneous Boltzmann equation by 
making the ansatz /i(x, k, t) = A^B(ek-T„vi(t)/n + V{:x. - 
xi(t))), with a similar expression for /o(x, k, i). This 
ansatz is now parameterized by the centcr-of-mass ve- 
locities vo.i(i) and positions Xo,i(i) of the two atomic 
clouds. From this, we get the equations of motion. 



iVm— ^ = 2Nmn x vi 
dt 
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Nm^^ = 2Nmn x vq 
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(6) 



with A'^ the particle number per spin state. Note that due 
to the centrifugal force, we need to have that \fi\ < ujq. 

We again linearize the above equations using that 
r(v, x) ~ r'v due to the isotropy of the collision in- 
tegral. We next observe that all the motion in the z- 
direction decouples. We therefore only consider the mo- 
tion of the clouds in the x — y-plane, since this contains 
the spin drag Hall effect. The linearized equations then 
yield a collective-mode spectrum with eight modes in to- 
tal. There are four modes in which the two clouds of par- 
ticles move in phase, and in which there is, as a result, no 
drag effect. The modes correspond physically to in-phase 
harmonic oscillations of the two clouds with frequencies 
wq ± SI. There are two different frequencies because the 
degeneracy due to the two equivalent directions of os- 
cillation in the effective two-dimensional system, is split 
by the external rotation. The four out-of-phase modes 
correspond physically to the two atomic clouds moving 
relative to each other. This results in transfer of momen- 
tum between the two clouds, leading to spin drag and 
damping of these modes. These modes have the frequen- 



cies 



to = — Z7 ± SI - 



{ij ± ny 



(7) 



The imaginary part of the above frequencies gives the 
damping rate of the modes, and is for O <C wg given 
by 7 = 1/2t = V /Nm with r the spin drag relaxation 
time. This relaxation time gives the longitudinal spin 
drag resistivity as pxx = Pyy = m/nr, and is estimated 
next. From the eigenvectors of the modes we find that the 
plane of oscillation of the out-of-phase dipole mode is not 



fixed in the co-rotating frame, which implies a transverse 
spin current. This is the trap equivalent of the spin drag 
Hall effect discussed in the previous section. 

Spin drag relaxation time — In the inhomogeneous 
case, we find that 



hk 



r(vo-vi,xo-xi)=ydxy — 3, 

NBiek-m^,{t)/h + V{x-MtM , (8) 



where 

rcon[./o,/i] = — ^(^01 )^ 



dko 



dk. 



dkA 



h '^"' ' J {2tt)3 J {2tt)3 J (27r)3 
xS{k + ks - kg - k4)5(ek + ek2 - ^ka - <^kj 
x{[l + /i(x,k,i)][l + /o(x,k2,t)]/i(x,k3,i)/o(x,k4,i) 

-/l(x,k,t)/o(x,k2,i)[l + /i(x,k3,i)][l + /o(x,k4,i)]}. 

Here, Tq^ is the two-body T-matrix, which equals 
ATTafi^/m, with a the scattering length for intcr-spin- 
species collisions. Introducing the response function 

X(x;q,w) (9) 

dk iVs(ek + l^(x))-7VB(ek+q + y(x)) 



(27r)3 



we find 
F' 



2B\2 



-xnr) 



Ek+q - e^ + TiLO + i{) 
dq 



Unfi' ""' ' J ""J (27r)3 
[Im[x(x;q,tj)]]^ 
sinh^ {l3huj/2) ' 



dug 



(10) 



from which we can determine the spin drag relaxation 
time. The imaginary part of the response function is 
worked out explicitly to yield 



Im[x(x;q,w)] 



log 
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with A the thermal de Broglie wavelength and /i(x) = 
t^-V{x). 

We estimate the above expression for a three- 
dimensional homogeneous system with density n for 
which we, in first approximation, have take the central 
density in the trap to make connection with the inho- 
mogeneous case. The result for r is evaluated numer- 
ically, and is shown in Fig. [T] We see that Bose en- 
hancement is indeed at play: 1/r gets dramatically bigger 
as the temperature approaches the critical temperature. 
We do find, however, that it remains finite as T | Tc- 
Furthermore, from our numerical results, we find that 
1/t(T) - l/r(r,) oc - V^^ oc T, - T. 
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^. Considering that the trapping 
0.01 — 1 kHz, this means the 



FIG. 1: Bose enhancement of the spin drag relaxation time 
r upon approaching the critical temperature T^ for Bose- 
Einstein condensation from above. 



Discussion and conclusion — We have introduced the 
spin drag Hah effect, i.e., the emergence of a transverse 
spin current, in rotating Bose mixtures. To determine 
whether the spin drag relaxation rate can be measured 
in principle, we make estimates of the relaxation time 
for some realistic values of the parameters. Taking for 
instance ^^Rb at temperatures between 10 and 100 nK, 
with an inter-species scattering length of about 100 Bohr 
radii, we find values of the order of 1 — 100 ms for densi- 
ties n = 10" -10^2 cm 
potential usually has ujq 
damping should indeed happen on an observable time 
scale. To compare with electronic systems we note that 
the Drude formula ni/ne^T (with e the electronic charge 
to convert to units of electrical resistivity) with our re- 
sult for r yields resistivities of the order of 10~^ — 10~^ 
r2m, many orders of magnitude larger than the spin drag 
resistivity in an electronic system |17| . 

One interesting aspect of the Bose enhancement of the 
1/r is the behavior close to the critical temperature. 
Here, we numerically found that 1/t(T) — l/r(Tc) ex 
— (T—Tc)'^, with K = 1 within our numerical accuracy. In 
future work we intend to investigate the value of this ex- 
ponent in various dimensions with renormalization-group 
methods, taking into account critical fluctuations that 
are not captured by the Boltzmann approach presented 
here. Note that in our previous work concerning one spa- 
tial dimension, we found a di verg ence of the spin drag 
relaxation rate as {T~Tc)~^^^ [13[, also with Boltzmann 
methods. An interesting aspect of a two-component Bose 
mixture is that it also may become ferromagnetic above 
the critical temperature of Bose-Einstein condensation. 
We intend to study also the effects of this transition on 
the spin drag. 

The collective mode spectrum determined theoretically 
in this Letter can be observed experimentally by setting 
the two spin states in relative motion. This can for exam- 



ple be achieved by shortly applying small magnetic field 
gradient to excite the spin-dipole mode. Another possi- 
bility is to use a state-selective laser in a manner that 
is similar to the generation of the second-sound dip ole 
mode in a partially Bose-Einstein condensed gas [18 1. 



We hope that the close collaboration between theory 
and experiments in this area, will lead to more insight 
into bosonic transport and, on the long run, may eventu- 
ally lead to the development of useful atomtronics devices 
[l9| , where atoms rather than electrons are the main car- 
riers of transport. 
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